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Replica symmetry breaking solution for the fermionic Ising spin glass and the 

Ghatak— Sherrington model 

H. Feldmann and R. Oppermann 





We solve the fermionic version of the Ising spin glass for arbitrary filling fi and temperature 
T taking into account replica symmetry breaking. Using a simple exact mapping from fi to the 
anisotropy parameter D, we also obtain the solution of the 5=1 Sherrington-Kirkpatrick model. 
An analytic expression for T = gives an improved critical value for the first-order phase transition. 
7j , We revisit the question of stability against replica-diagonal fluctuations and find that the appearance 

/^S ■ of complex eigenvalues of the Almeida-Thouless matrix is not an artifact of the replica-symmetric 



approximation. 



I. INTRODUCTION AND MAPPING OF THE TWO MODELS 



A direct generalization of the Sherrington-Kirkpatrick (SK) model M to include quantum fluctuations are fermionic 
spin glasses g]. They provide a larger class of models, since the partition functions of all classical spin glass models 
can be extracted from them through the Popov-Fedotov trick [00. They also allow to extend standard spin glass 
theory to itinerant systems, to study the influence of spin glass order on the excitation spectrum [^j and to investigate 

C/3 I the competition of spin glass order with other kinds of ordering typical of quantum systems Q . All these aspects may 
be relevant to the physics of heavy fermion compounds. 

Another extension of the SK model is the S — I spin glass in a crystal field, realized for example in (Tii_xVx)203 
M. Both extensions show tricritical behavior as the chemical potential fi (for fermionic spin glasses) or the anisotropy 
parameter D are varied [[7|,y|. A simple mapping relates the two models, as far as the static properties are concerned. 

fH , The fermionic Ising spin glass (ISG/) is described by the grand canonical Hamiltonian 

o; 

Oj! H = -'^Jijaiaj - fj.'^n^ (1) 



^ . The coupling Jy is Gaussian distributed around with variance J^. Below, we always set J — 1. The main difference 
TJ" ■ to the SK spin glass is that spins and occupation numbers are given in terms of fermionic operators which act on a 
Cn , space with four states per site, |00), | t 0), |0 j), and | tj): 

O . 

^•^ I a — ala^ — ajax n ~ ala^ + ajox (2) 

^\^ To obtain the thermodynamic behavior of the model, we calculate the free energy using the replica trick |9[| . Integrating 
"^ I over the distribution of Jij creates eight-fermion correlations, which are decoupled using the Parisi matrix of order 
C^ I parameters Q"'^''''^ . For mean field theory, we use the static saddle point of this matrix. Details of this calculation 
can be found in |g,|l^ . 

The S — 1 anisotropic spin glass or Ghatak-Sherrington (GS) model JTl is represented by the Hamiltonian 



\ ^ 7 c c in \ ^ c2 



O : ^^-l^J^i S^.S,. +D}_^St (3) 

.5^ ' where 5*2 may have the values —1,0, and 1. It is easy to see that the partition function corresponding to the 

^>^ . Hamiltonian (|l|) is identical — apart from the constant 1 + exp(— 2/3/j,) — to the one deflned by Eq. (g), provided one 

^H maps anisotropy and chemical potential according to 



aPD _ ,/3m j_ ,-l3t^ (4) 



e'^ — e '^ + e 



As a consequence of (y), the thermodynamic properties of both models are directly related. At T — > 0, one may even 
neglect the last term and set Z? = /x. Note, however, that there is a whole class of fermionic correlations that can 
not be expressed in terms of Sz and are therefore unique to the ISG/, because their definition requires a fermionic 
generating functional. A good example is the fermion Green's function llj]. Also, the Heisenberg and XY variant of 
the ISG/, since they exhibit quantum dynamics, do not have a direct classical analogue. 



In passing, we note that there is another interesting connection between the two models in addition to the mapping 
(y) and the above mentioned Popov-Fedotov method. One could introduce a repulsive Hubbard-like on-site interac- 
tion U to Hamiltonian (||), which would deplete the doubly occupied states in the limit of infinite strength. In this 
case, the Hamiltonian and a mapping similar to (m would read 



n = -^Jijdiaj - {^i-—)^n^ + U^{a\a^ - -){a 



lax 



and 



e^^ = e 



-Pf, 



„/3m-/3(7 



(5) 



II. RESULTS AT NONZERO TEMPERATURE 

It has been known for a long time that the two models under consideration are described by two mean-field 
parameters, q{x) and the replica-diagonal saddle-point value, Q"" = q (in the notation of the ISG/). 

Already from replica-symmetric (RS) calculations, one can locate the second-order critical line, terminating at the 
tricritical point at (T = 1/3, D = 0.962 or /.i = 0.961). In order to correctly describe the system below this line, 
one has to allow for replica symmetry breaking (RSB). Since only few results are available for full RSB |6), the first 
step of the Parisi scheme (IRSB) [p2l can give valuable information. In many cases, it already comes close to the full 
solution or allows a good conjecture. 

Using the standard form of RSB at the one-step level, we use the order parameters qi = q{l) and 52 = 9(0). The 
position of the jump in q is given by to. After integration over the fermionic fields and some other transformations 
we obtain the replica-broken free energy 



f3f =^((9 - 1)' - (qi 1)' + m{q! <il)) - /3m - - 

4 m 

C = cosh{l3h + pj^qi - q2Zi + (3J^Z2) + cosh(/3//)e-^ («--?!) 



G j-G 

In / C™ with 



(6) 

(7) 



The mean-field solutions of the model are given by the condition that / is stationary with respect to g, qi, (72, 
and TO,. Considerable numerical effort is necessary to solve the resulting self-consistency equations for these four 
parameters simultaneously, and for all values of chemical potential and temperature within the ordered phase. In 
Fig.|l| we present a selection of our results, introducing the susceptibility x = (3{q — (1 — m)qi — mq^). We compare 
our results with results from the replica-symmetric calculation. Note that x in IRSB is already very close to the 
expected exact value x = 1 in the ordered phase, especially for temperatures close to the phase transition. 
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FIG. 1. Susceptibility x fi'Hd Edwards-Anderson order parameter q(\) for the ISG/ for half-filling and for \i — 0.5. Solid and 
dashed lines show the results from the 1-step replica broken and the replica-symmetric solutions, respectively. 



III. ZERO TEMPERATURE RESULTS 



In the limit of vanishing temperature, q = q + 0(T). Also to vanishes as T — > 0. Therefore, for the zero-temperature 
behavior of the model, we have to use a new set of parameters in order to avoid divergences in Eq. (|6|) . We replace 
f}J{q — qi) by the single-valley susceptibility x and m by aT. Note that although the to, the size of the "lower step" 
in q{x), vanishes with the temperature, replica symmetry breaking still has a profound effect on the system. 



In these new variables, the Umit of Eq. (^) can be written as 

/ = Uxiqi - 1) + ^Mqj -ql)-f^-- f In/ 
where / is given for the case ^/J < x/2 by 

2 ^ \/gT 



(8) 



_|_ _g-a^q^Z2 + ^a'^{qi-q2) fl _|_ gj.f 



. -^/g^Z2 + Q(gl - g2) 

Vgi - g2\/2 



(9) 



and for the case /i/ J > x/2 by 



/ — _ga\/9^^2 + |a^(gi-g2)| 2^ _^ gj.f^ 



' (l + erf ( 



g2Z2 - (7 - I) + a(gi - g2) 



+ Ie''(7-I)(erf( 
1 



Vgi - g2\/2 



Vgi - g2V2 



)) 



g2Z2s 



Vgi - g2V2 



(10) 



ig-a\/92«2 + ^a^(gi-g2) (2^ 



erf 



■ -Vg2g2 - ( J - 2) + °(gl ~ 92) 
^gl - g2V2 



Again, variation with respect to the order parameters gives four equations describing the saddle point solution. We 
solved these for all chemical potentials within the ordered phase and present the results for q = qi, x and x in Fig. 0, 
together with the replica symmetric order parameters for comparison. In RS approximation, ^ = ;^, which illustrates 
the importance of RSB. 

We also calculated the free energy from Eqs. (H)-(|o|) and compared it with the free energy for the paramagnetic 
phase, /para — — 2/i, to obtain the first order phase transition at /it,i — 0.881, a slightly lower value compared with 
the replica-symmetric approximation, which gave /it,o — 0.900. 

The transition between the two regimes described by Eqs. {m and {Km takes place at fi — 0.119, which is also 
related to the width of the band gap at half filling |l^,|l^ . 
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FIG. 2. Normal (x) and single-valley (x) susceptibility (left) and spin-glass order parameters (right) at zero temperature. 
Solid curves correspond to the replica broken solution, dashed curves to the replica-symmetric one. For each case, the left 
vertical line indicates the onset of the //-independent solution for small fi, while the right vertical line locates the first order 
transition to the paramagnetic regime. 



IV. REPLICA SYMMETRY BREAKING IN CHARGE CORRELATIONS 



As sections O and p^ show, RSB has a profound effect on the low temperature phase. However, spin and charge 
correlations are affected differently. While spin correlations are modified by RSB within the whole ordered phase (see 
the susceptibility in Figs, n^ and 0), charge correlations remain essentially unchanged by RSB within part of the phase 
diagram. 



The most obvious observable parameter connected to the rephca-diagonal order parameter q is the filhng v. Both 
are related to all orders of RSB by 

q = l-coih{(3^i){v-l) (11) 

Zero-temperature results for v are displayed in Fig. g. The difference between the broken and the unbroken solution 
is largest at the point where the RS value of u turns constant. Around /i = 0.7, the IRSB solution bends around to 
come close to the RS one and the two solutions cross. But the difference stays comparatively small, so for large jjl 
the filling is virtually not affected by RSB. A similar situation appears in the standard spin glass theory for to > 1 
component spin glasses. There, in a strong magnetic field two characteristic lines appear, first the Gabay- Toulouse 
line M-na], where RSB appears for the order parameter perpendicular to the field while it has little effect on the 
longitudinal order parameter. At a lower temperature, close to the line found by de Almeida and Thouless, also the 
longitudinal order parameter strongly feels RSB. In this context, one may view the charge and spin degrees of freedom 
to be "orthogonal" and the chemical potential to correspond to a magnetic field. 
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FIG. 3. pi - dependence of the filling at zero temperature in RS and IRSB approximation. For fj, < 0.119, both approximations 
yield a constant u, while the full solution, estimated with the dashed line, gives u — 1 only at /i = 0. Note the remarkable drop 
of the RSB solution to the RS line around fi = 0.7. After the crossing the two lines stay close together. The inset shows the 
difference between the IRSB and the RS solutions for large /i. 



For nonzero temperatures, the difference between RS and IRSB in general is smaller than for T = 0, as shown in 
Fig. H. The merging of the two lines becomes smoother with increasing temperature. The effect of broken replica 
symmetry on i^ already is invisible at T = 0.4 in Fig. U, far below the second order phase transition (see Fig. ^. 
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FIG. 4. Filling factor u as a function of the chemical potential for temperatures T — 0,0.05, 0.1 ... 0.4. For each pair of 
curves, the lower one gives the RS solution, the upper one the IRSB result. We applied an offset in steps of 0.1 to separate 
the results for different T. For fj, — 0, there is always half-filling. Note that the effect of RSB on i^(^) is significant only in a 
intermediate range of chemical potentials and for small temperatures. 

For the GS model, the results of this section on the filling can be directly translated to the average number of sites 
with Sy = 0: 



(s! 



a -Si 



coth{Pfi){iy - 1) 



(12) 



V. REPLICA-DIAGONAL STABILITY 



The replica-symmetric solution is unstable against replica-symmetry breaking, as shown by de Almeida and Thou- 
less (AT) in [|l6| . This problem, together with its famous solution using an ultrametric saddle-point matrix ||9|, is one 
of the main reasons for the enormous theoretical interest in spin glasses over the past decades. 

In addition to the eigenvalue that marks the onset of replica-symmetry breaking, AT obtained two additional pairs 
of eigenvalues which merge in the replica limit. For the Sherrington-Kirkpatrick model, these were positive for all 
temperatures and did not pose a problem. In the case of the Ghatak-Sherrington model and the ISG/, however, the 
replica-diagonal eigenvalues can become complex, as noted by several authors p7[- p^ , p^ . 

This result is very difficult to interpret. Lage and de Almeida |1^ derived another set of stability conditions and 
found their condition d^f/dq^ to be violated at low temperatures, in a region determined by the "crossover line" of 
Ref. Ig]. Mottishaw and Sherrington |1^ pointed out that the system is unstable against RSB and they suspected 
that the full Parisi solution will have only real eigenvalues again. 

We have numerically evaluated the self-consistency equations from Ref. B within the ordered phase and calculated 
the replica-diagonal eigenvalues. Directly below the second-order phase transition and at low temperatures, they are 
real and positive, indicating stability against replica-diagonal fluctuations. But in between there is a region extending 
from fi — to the tricritical point, where they have an imaginary part, while the real part remains positive (see Fig. 
ph. At half-filling, the complex eigenvalues occur between T — 0.53 and T = 0.61. 
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FIG. 5. Region of complex replica-diagonal AT eigenvalues in the replica-symmetric approximation. Also shown is the line 
of continuous transitions ending in the tricritcal point. 



To extend the AT analysis of the replica-diagonal eigenvalues to replica symmetry breaking, we start with the 
ansatz for the eigenvector 



|g(aa)| 
|^(afc)| 



(13) 



where the replica indices run from a,b = 1 to n and only pairs with a < b are considered for 77^°''). For k step RSB, 
the replica-diagonal eigenvector of AT can be generalized to 



:(''") = a and t]^"''^ 



Pi if Q"-'' = qi 

(32 if Q^" = 92 



(14) 



Substituting Eq. dlj) into the quadratic form 



^i^G^l with G(ab)M) ^ J , J , / 



dQ"^ dQ 



(15) 



yields fc + 2 replica-diagonal eigenvalues. In the limit of n — ^ 0, these eigenvalues can also be represented as the 
eigenvalues of a matrix of (much smaller!) size k + 2: 



Mk = 



f" 



/'i 



-2 /•/ -2 f 

(l-mi)-' 1 (1-mi)-'- 



11 



(mi— rM2) 



/'2 



(mi— 7712) 



-2 



.,)f'k 



(m(fc-l)-™k) 
1 -2 f, 

Here, we have introduced the shorthand notations 



21 



r/fei 






f'k+l 



-2 f 

(mi— m2) .'2(A;+1) 



-2 



r/fc 



(m(fc_i)-mfc) ^'^(fc+l) 
mfc 



(16) 



-/(fe+i)(fc+i) / 



d'f 



d df 



d df 



f" — f' = — and f = ■' 

9(7^ ' ' dqi dq ' *"' 9(7i Sg^ 



(17) 



We have carried out the above calculations numerically for IRSB with the order parameters presented in section O. 
Regions with complex replica-diagonal eigenvalues are shown in Fig. g. Now there appear two of them. Since / and its 



derivatives are real, obviously complex eigenvalues always have to appear in pairs conjugate to each other. Labelling 
the three eigenvalues at IRSB in proper order, we find that in the high-temperature region of Fig. ^ Ai and A2 turn 
complex, while in the low-temperature region A2 and A3 have nonzero imaginary parts. Therefore, the two regions 
can be clearly distinguished. The high-temperature region has approximately the same shape as the one from RS 
calculations. Its boundary at /i = is given by T = 0.52 and T — 0.61. This fact is a strong indication that the effect 
of RSB is small in the corresponding temperature range and that the occurance of complex AT eigenvalues continues 
for arbitrary RSB. Complex eigenvalues thus seem to be a consequence of the replica limit. Their interpretation as 
far as stability is concerned remains an open question. 




FIG. 6. The gray regions show complex AT eigenvalues in one-step replica symmetry breaking. Dashed boundaries indicate 
difficulties in the numerical algorithm. Here, the true boundary may be slightly shifted. 

VI. CONCLUSION 

Using the equivalence of the thermodynamic properties of the Ghatak-Sherrington model and the fermionic quantum 
Ising spin glass, we simultaneously solved these two models in the first step of the replica symmetric approximation 
combining analytical and numerical methods. We also obtained the T = limit of these solutions. As expected, the 
effects of RSB are comparatively large in general, while the order parameters of charge are almost unchanged by RSB 
in certain regions of the ordered phase. Analyzing the replica-diagonal AT eigenvalues, we found that the puzzling 
complex stability eigenvalues, previously found in RS calculations, also appear in one step RSB. They even cover a 
larger part of the phase diagram. On this basis, we formulated the conjecture that complex eigenvalues of the stability 
matrix appear at any order of RSB. 
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[1] D. Sherrington and S. Kirkpatrick, Phys. Rev. Lett. 35, 1792 (1975). 

[2] R. Oppermann and A. Miiller-Groehng, Nuc. Phys. B 401, 507 (1993). 

[3] V. N. Popov and S. A. Fedotov, Zh. Eksp. Teor. Fiz. 94, 183 (1988), [Sov. Phys. JETP 67, 535 (1988)]. 

[4] O. Veits, R. Oppermann, M. Binderberger, and J. Stein, J. Phys. I France 4, 493 (1994). 

[5] R. Oppermann and B. Rosenow, cond-mat/9803249 (1998). 

[6] H. Feldmann and R. Oppermann, Eur. Phys. J. B 10, 429 (1999). 

[7] S. K. Ghatak and D. Sherrington, J. Phys. C 10, 3149 (1977). 

[8] B. Rosenow and R. Oppermann, Phys. Rev. Lett. 77, 1608 (1996). 

[9] M. Mezard, G. Parisi, and M. A. Virasoro, Spm glass theory and beyond (World Scientific, Singapore, 1987). 



[10] 


R. 


[11] 


R. 


[12] 


G. 


[13] 


M, 


[14] 


M, 


[15] 


D. 


[16] 


J. 


[17] 


E. 


[18] 


P. 


[19] 


F. 



Oppermann and B. Rosenow, to appear in Phys. Rev B (1998). 
Oppermann and B. Rosenow, Phys. Rev. Lett. 80, 4767 (1998). 
Parisi, J. Phys. A 13, 1101 (1980). 

. Gabay and G. Toulouse, Phys. Rev. Lett. 47, 201 (1981). 
. A. Moore and A. J. Bray, J. Phys. C 15, L301 (1982). 
M. Cragg, D. Sherrington, and M. Gabay, Phys. Rev. Lett. 49, 158 (1982). 
de Almeida and D. Thouless, J. Phys. A 11, 983 (1978). 

J. S. Lage and J. R. L. de Almeida, J. Phys. C: SoUd State Phys. 15, L1187 (1982). 
Mottishaw and D. Sherrington, J. Phys. C: Sohd State Phys. 18, 5201 (1985). 
A. da Costa, C. S. O. Yokoi, and S. R. A. Salinas, J. Phys. A 27, 3365 (1994). 



